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Abstract 

The screened QED corrections of the first orders in a and 1/Z to the g factor and the hyperfine splitting 
of lithiumlike ions are evaluated within ab initio quantum electrodynamical approach. The complete gauge- 
invariant set of the two-electron self-energy diagrams in the presence of the magnetic field and a dominant 
part of the two-electron vacuum-polarization diagrams are calculated. The most accurate values of the g 
factor of Li-like lead and uranium are presented. The theoretical prediction for the specific difference of the 
hyperfine splittings of H- and Li-like bismuth is improved. 
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I. INTRODUCTION 



The bound-electron g factor has been the subject of intense experimental and theoretical inves- 
tigations over the past decade. Recent measurements for \ow-Z hydrogenlike ions with a spinless 
nucleus have reached the precision of 1CT 9 [1-3]. Together with the corresponding theoretical 
studies these experiments have lead to the new value of the electron mass, four times more ac- 
curate than the previously accepted value (see Ref. [4] and references therein). Experimental 
investigations of ions with more than one electron are anticipated in the nearest future. In particu- 
lar, measurements of the g factor of H-like and Li-like calcium and silicon are currently in progress 
by the Mainz-GSI collaboration [5]. An extension of these studies to high-Z H-like, Li-like, and 
B-like systems planned in the framework of the HITRAP project [6, 7] will provide a stringent 
test of the bound-state QED in the strong electric field of the nucleus. Moreover, investigations 
of the g factor of heavy B-like ions can lead to an independent determination of the fine structure 
constant [8]. The motivation for studying the g factor of Li-like and B-like ions follows from the 
higher theoretical accuracy that can be reached for a specific difference of the g factor values of 
H-like and Li-like ions (or H-like and B-like ions) of the same isotope. Various effects on the 
g factor of H-like ions were investigated during the last two decades: one-loop [9-16] and two- 
loop [17, 18] QED corrections, recoil corrections [19-21], nuclear polarization effect [22], and 
nuclear size effect [23]. The theoretical investigations of the g factor of Li-like ions were con- 
ducted in Refs. [24-27]. Apart from the one-electron contributions to the g factor of 2s state, the 
effects of the interelectronic interaction should be taken into account in three-electron ions. The 
one-photon exchange correction was evaluated in the framework of QED in Ref. [25]. The higher- 
order contributions of the interelectronic interaction were calculated by means of the large-scale 
configuration-interaction Dirac-Fock-Sturm method in Ref. [26]. Still, the uncertainty associated 
with these contributions amounts to more than half of the total theoretical uncertainty. The ef- 
fect of the interelectronic interaction on the QED corrections was treated within two approaches. 
For low-Z ions the perturbation theory to the leading orders in aZ was employed [24, 26]. For 
middle-Z and high-Z ions more accurate results were obtained by evaluating the one-electron 
QED corrections in an effective screening potential [27]. Nevertheless, for all values of Z the 
uncertainty of the screened QED effects contributes significantly to the total uncertainty, and the 
rigorous evaluation of these effects is in demand. 

Hyperfine structure of highly charged ions comprises another sensitive tool for probing QED 
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effects in strong fields. Accurate measurements of the ground-state hyperfine splittings were per- 
formed for several H-like ions, including 209 Bi, 165 Ho, 185 Re, 187 Re, 207 Pb, 203 T1, and 205 T1 in Refs. 
[28-32]. These experiments motivated corresponding theoretical investigations [33-40]. The the- 
oretical uncertainty of the hyperfine splitting is also dominated by the nuclear effects, mainly by 
the Bohr-Weisskopf effect [35]. The theoretical investigations have shown [40] that simultaneous 
studies of the hyperfine splitting in H-like and Li-like ions of the same isotope can significantly 
improve the accuracy of the theoretical prediction. As a result, the ground-state hyperfine splitting 
in Li-like bismuth was predicted to a high accuracy using the experimental result for the Is hy- 
perfine splitting in H-like bismuth [38-40]. The indirect measurement of the hyperfine splitting of 
lithiumlike bismuth performed in Livermore yielded value of 820(26) meV [41]. Determination 
of this splitting to a much higher accuracy (of about 1CT 7 ) is planned at GSI in the framework of 
the HITRAP project [42]. This requires further improvements of the theoretical predictions for 
Li-like ions and, in particular, evaluations of the QED screening effect. An approximate treatment 
of this effect was accomplished in Refs. [39, 43-47] by employing an effective screening potential 
in calculations of the one-loop self-energy and vacuum-polarization diagrams. 

Rigorous evaluation of the two-electron self-energy and vacuum-polarization diagrams (Figs. 
1 and 2) remained a challenge for theory until recently. In our Letter [48] the complete aZ- 
dependent contributions of the two-electron self-energy diagrams and a dominant part of the two- 
electron vacuum-polarization diagrams have been evaluated for the hyperfine structure of Li-like 
bismuth and for the g factor of Li-like lead. In the present paper we describe in detail the evalua- 
tion of the screened quantum electrodynamical correction in presence of external magnetic field. 
Furthermore, we extend our calculations to the wide range of nuclear charge Z = 20 - 83 in case 
of the hyperfine splitting. The accuracy of the theoretical prediction for the specific difference of 
the hyperfine splittings of H- and Li-like bismuth is improved. As to the g factor, we present the 
results for lead and uranium ions. 

The relativistic units {h—\,c — \,m — Y) and the Heaviside charge unit [a = e 2 / (47r), e < 0] 
are used throughout the paper. 



3 



II. FORMULATION 



A systematic derivation of the QED corrections in a fully relativistic approach requires the use 
of perturbation theory starting with a one-electron approximation, described by the Dirac equation, 

[-ia • V + /3m + V(x)] ^ n (x) = e re ^„(x) . (1) 

In our present treatment the binding potential V(x) = V(|x|) denotes the nuclear potential only. 
The interaction of the electrons with the quantized electromagnetic field and the interelectronic- 
interaction effects are accounted for by the perturbation theory. In this way we obtain quantum 
electrodynamics in the Furry picture. To derive the formal expressions for the perturbation theory 
terms, we employ the two-time Green-function method [49]. 

The diagrams, which contribute to the screened self-energy and vacuum-polarization correc- 
tions of the first order in a and 1/Z in presence of external magnetic field are depicted in Figs. 
1 and 2. In order to simplify the derivation we specify the formalism, where the electrons of the 
closed shell are regarded as belonging to a new redefined vacuum. The redefinition of the vacuum 
results in replacing iO by — iO in the electron propagator denominators corresponding to the closed 
shell. In this formalism the one-electron radiative corrections are incorporated together with the 
interelectronic-interaction contributions. In particular, the one-loop two-electron contributions are 
merged with the two-loop one-electron contributions. The corresponding diagrams are depicted in 
Fig. 3. Below we briefly describe the scheme of derivation of the formulas, corresponding to this 
set of diagrams within the two-time Green-function method. In order to obtain the two-electron 
corrections one may simply consider the related expressions with the standard definition of the 
vacuum and then consequently make a replacement 

J2^^2,i5(e-e c )J2\c)(c\ (2) 

n n c 

for each of the electron propagators inside the loops. Here and in the following the notation 
u = 1 — iO is used. The summation over c is performed over all electrons of the closed shell. 

To zeroth-order approximation the state \a) of the electron is defined by the Dirac equation (1). 
The energy shift of an isolated level due to the interaction is given by [49] 

(2m)-^ r dE(E-Ej 0) )Ag aa (E) 

In our case the unperturbed energy E^ is the Dirac energy e a from Eq. (1). The contour Y 
surrounds only the pole E = E^ , Ag aa (E) = g aa (E) - g { aJ (E), g aa (E) = (ip a \g(E)\tp a ), and ip a 



is the unperturbed wave function. The time-Fourier transform of the Green function is defined as 

i r+oo 

g(E, x', x)5(E — E') = / dtdt' exp(iE't' - iEt)(0\Tip(t' ,x'UUt : x)\0) . (4) 

The Feynman rules for the Green function are given in Ref. [49]. The energy shift AE a and the 
Green function g(E) are to be expanded as a power series in a, 

AE a = AEW + AE® + AE® + . . . , 
Ag aa (E) = AgU(E) + Ag^(E) + Ag^(E) + .... 

Then from Eq. (3) we find the first-, second- and third-order terms of the energy shift, 
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Equation (9) and the Feynman rules for g(E) yield the formal expressions for the total contribution 
of the two-loop diagrams presented in Fig. 3. Consequently replacing each electron propagator 
according to Eq. (2) we obtain the formal expressions for the two-electron one-loop diagrams 
displayed in Figs. 1 and 2. 

For brevity we introduce the operator 

I(u) = e 2 a»a u D^{uj) , (10) 



where a M = (1, a) are the Dirac matrices, and is the photon propagator. It is given by 

exp(iu)|x 12 |) 



D^u^xia) = g^- 



47r|xi 2 | 



(11) 



in the Feynman gauge and by 

A)o(w,Xi 2 ) 



47r|xi 2 | ' 



Ao(w,xi 2 ) = £> fc(w,xi 2 ) = 0, 



n ( x r exp(^|x 12 |) 1 - exp(iw|xi 2 | 

Afc(w,x 12 ) = d ik — | h ViiV 2fc - 



AlT X12 



47ra; 2 |xi 2 | 



(i, k = 1,2, 3), (12) 



in the Coulomb gauge. Here x 12 = x x — x 2 , Co = y/cu 2 + iO. The branch of the square root is 
fixed by the condition Imw > 0. In order to handle the infrared divergencies it is convenient to 
introduce the photon mass fi. In the Feynman gauge, it results in the replacement \/uj 2 + iO — > 



The interaction with the external magnetic field can be represented by an operator T . For both 
cases under consideration (g factor and hyperfine splitting) it is proportional to [r x a] 2 , what 
defines the angular momentum structure. Explicit formulas for T will be given in Sec. V. 

III. SCREENED SELF-ENERGY 

The diagrams of the screened self-energy correction, corresponding to the first term in Eq. 
(9), are shown in Fig. 1. We divide the total contribution of these diagrams into the reducible 
and irreducible parts. The irreducible part is the sum of the terms where the energies of the 
intermediate states are different from the energy of the initial state. The reducible part is the 
remainder. We denote the irreducible parts of each diagram A - F by the same letter: A£gQ^~ F \ 
It is convenient to divide them into three groups, according to the number of the ^-dependent 
denominators (u> is the virtual photon energy, over which the integration is performed). The terms 
with only one denominator (A, B and E) are referred to as a "modified self-energy" terms, since 
all of them have the form of a matrix element of the self-energy operator, (X\T,(e) \Y). The terms 
with two denominators (C and F) are denoted as a "modified vertex", and the diagram D with 
three denominators is denoted as a "double-vertex". The reducible parts are considered together 
with the related contributions that arise from the second and the third terms in Eq. (9). Their sum 
is divided into three parts: G, H, and /, according to the number of the ^-dependent denominators. 
These three parts are considered together with the "modified self-energy", "modified vertex", and 
"double-vertex" terms, respectively. Finally, the total contribution of the two-electron self-energy 



a/w 2 — + iO. The limit /i — > should be taken after 
I(cu) has the following symmetry properties: 



removing the divergencies. The operator 



I(u>) = I(-cu) , 




(13) 
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diagrams is given by 

= £ EM)"* H«Td> + Alggg + A£-g + Aifgg 

i A pSE(_B) a pSE(F) a pSE(G) . a pSE(H) » pSE(J) \ 
+ Z ^- fc SQED + ^-^SQED + ^-^SQED + ^-^SQED + ^-^SQED J ■ U^J 

Here P and Q are the permutation operators, interchanging the valence (a) and the core (b) electron 
states, (— l) p is the sign of the permutation P. The summation over b runs over two core electron 
states with different projections of the total angular momentum. In what follows, we will also use 
the notation A = EQb — Spb- 

One has to pay special attention to the ultraviolet (UV) divergencies, that arise in the formal 
expressions under consideration. First, we introduce the unrenormalized self-energy operator, 

n 



Here and below the integration over u is carried out from — oo to +00. Every diagram involving 
the self-energy loop has to be considered together with the corresponding diagram with the mass 
counterterm, what results in the replacement E(e) — > S R (e) = ^°5m. The matrix elements 

of £r(£) still have the divergent part: 

(p\E R (e)\q) = [e - a ■ p - f3m - V(x)] \q) + finite part , (16) 

where B^> is the UV-divergent constant. Assuming that \p) and \q) obey the Dirac equation (1) 
we have, 

(p|S R (e)|g) = B {l \e - e p )5 pq + finite part. (17) 

In order to isolate the divergent part we follow the potential-expansion approach [50]. The fi- 
nite part of the self-energy matrix element is then divided into the zero- and one-potential terms 
evaluated in momentum space, and the many-potential term evaluated in coordinate space. The 
diagrams with one vertex inside the self-energy loop also suffer from UV-divergencies. It can be 
shown that for an arbitrary operator U, 

— [duY , W^MM)^!^) LWm ) + m (18) 

2-k J ^ (e -u - ue ni )(e -u - ue n2 ) 



"1,2 v 



In our case U is either Tg or /(A) (in the latter case (p|{7|(7) = (pr|/(A)|q , s)). Due to the Ward 
identity we have = —B^\ The finite part of the vertex contributions is divided into the zero- 
potential term evaluated in momentum space and the many-potential term evaluated in coordinate 
space. 



From Eq. (17) one can see that for the first-order self-energy correction the divergent part is 
zero. However, in case of the higher-order diagrams under consideration the contributions of par- 
ticular diagrams are divergent. Nevertheless, as it is shown below, the sum of all the contributions 
to the screened self-energy correction is finite. 



A. "Modified self-energy" diagrams 

The irreducible parts of the diagrams A, B and E can be presented as the matrix elements of 
the self-energy operator (X\Y,(s)\Y) with various wave functions (X\ and \ Y). The formulas for 
these contributions are as follows, 

A pSE(A) _ A pSE(Al) A F SE(A2) , . q , 

SQED - SQED + SQED > 

A pSE(Ai) V-'/P,IW, \i„ \ (ni\T \n 2 )(n 2 Pb\I(A)\QaQb) 

A ^SQED = 2 2^ ( Pa \WPa)\ni) 7- —T7- — v , (20) 

[£Pa — £ m ) \ £ Pa — £ n 2 ) 



"1,2 



A«» = 2^V a |E(^)| ni ) ( " 1 f 6|/(A)l "f >( " 2|r °' W , (21) 

{£Pa — £ni)\£Qa ~ £ n 2 ) 



"1,2 



^Pa £-ni &Pa £ n 2 



"1,2 



A tttSE(£;) _ a j-,SE(Bl) a j-,SE(E2) 

SQED — ^-^SQED + ^-^SQED i V Z -JJ 

A^ifc" =2^V a |E( £p „)|, il ) < T 0l " 2>< " 1 w' /(A " Q '; W - < 24 > 

1 ^Pa — ^ni ) [£Pb — £n 2 ) 



"1,2 



(£Pa — £"iJ(£Q6 — £ n 2 ) 



"1,2 

All the reducible terms of the similar structure are denoted as A_EgQ^§, 



A 77iSE(G) _ a rpSE(Gl) a rpSE(G2) « t-,SE(G3) 
A -^SQED - A -^SQED + A -^SQED + A -^SQED > 

A^sqed^ = -2^7 Z [<Po|S(ep a )|ni>(ni|ro|Pa>(PaP6|/(A)|QaQ6> 

rai l £ Pa £ "lJ L 

+ (Pa| £( £Pa ) |m) (mP6|/(A) |QaQ6) (Pa|T |Pa) 
+ (Pa|S(£ Pa )|Pa)(Pa|T |n 1 )(n 1 P6|/(A)|QaQ6)] , (27) 
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a ^sqed 2) =2^' Z kPal^epa^HmlTolPaHPaPbll'iA^QaQb) 

E-Pa &n\ - 



ni 



+{Pa\Z(e Pa )\n 1 ){n 1 Pb\I'{A)\QaQb)({Qb\T \Qb) - (Pb\T \Pb)J 



+2 (Pa\L{e Pa )\Pa) ^ 



(PaPbirjAjlrnQfyimlTolQa) 
{PaPb\I\A)\Qan l ){n l \T \Qb) 



£Qb — £ni 

r,SE(G3) 



(28) 



A£ S qed = (Pa\Z(e Pa )\Pa}(PaPb\I"(A)\QaQb}{(Qb\T \Qb} - (Pb\T \Pb)J . (29) 

Equations (19)-(29) possess ultraviolet (UV) divergences. Taking into account the mass countert- 
erm and employing Eq. (17) we find that AE^q^> nas a non-zero UV-divergent part, 

Aiggg(UV) = 2 B« ^WWWiWW . (30) 

By the end of the next subsection we will show that the sum of all the UV-divergent terms is zero. 
B. "Modified vertex" diagrams 

For the irreducible parts of the diagrams C and F we have, 

A 7T>SE(C) _ a t-iSE(CI) a rpSE(C2) ,~ 1 , 

SQED ^-^SQED + ^-^SQED > W l ) 

ae se(ci) = 2 j_ f duj £n 'iyr Pa (Pan 2 \I(oj) Kn 3 ) (m \T \n 2 ) (n 3 Pb\I(A)\QaQb} 

AjB se(C2) = 2— [ dJ n 'y a ( Pa ^|/(^)|n 1 n 3 )(n 1 P&|/(A)|n 2 g&) (n 3 |T |Qa) 
sqed 2yr y ( £pa _ w _ U£ni ) ( £Qa _ w _ M£n2 ) ( £Qa - £n3 ) ' 

' '•l ,2,3 

A£; se(f) = 2 ± f du £n y b (Pan 2 |/(a;)|n 1 ga)(n 1 P&|/(A)|n 2 n3) (n 3 |T |Qfe) 

Since these diagrams have one vertex inside the self-energy loop, the corresponding expressions 
have the following structure of the cu-dependent denominators: (Ai — w)~ 1 (A 2 — uo)^ 1 . All the 
reducible terms that have similar structure are denoted as AE^f^ , 

A^sqS = A^QED^ + A£ S S QED 2) + Aijgg" , (35) 
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27t7 ^ (e Pa -uj-ue ni )(ep a -uj-ue n2 y I v v , 

711,2 

f, ^ (Pan 2 |/(o;)|n 1 ga)(n 1 P6|/ / (A)|n 2 g6) 

SQED " 2n J ( £pa _ w _, ueni )( ego _ a; _, uen2 ) 

x((Q6|To|Q6>-(P6|r |P6>), (37) 
A^sqed 3) = 2^' ^ [(Pa|S / ( £Pa )|n 1 )(n 1 |T |Pa)(PaP6|/(A)|gag6) 

^Pa £n\ - 



+(Pa\T!{e Pa )\n 1 )(n 1 Pb\I{A)\QaQb){Pa\T Q \Pa) 



+ (Pa|S / (e Pa )|Pa)|2^' 



ni 



' (Pa|T |7i 1 )(7i 1 P6|/(A)|gaQ6) 

^Pa £-ni 

| (P6|r |n 1 )(Pan 1 |/(A)|QaQ6) ' 

+ (PaP6|/'(A)|gag6) ((Pa|T |Pa> + (Qb\T \Qb) - (Pb\T \Pb)^) J . (38) 

Equations (3 1)— (38) diverge both in ultraviolet and infrared regions. The UV-divergent terms are: 

A^S»(UV) = 2L»> E '(^,N)(»^|/(A)IQ a W ^ 

£Pa &ni 

ni 

Agg{g'(uv) = 2L m£ {PaPb]m]nm{niWa) . (40) 

ni £Qa £ni 

ap s s q e S(uv) ^^.r'MWMi, (41) 

APsqed'^UV) = L {1) (Pa\T \Pa)(PaPb\I'(A)\QaQb) , (42) 
APsqed 2) (UV) = L« (PaP6|/ / (A)|g a g6) (<Q6|T |Q6> - (Pb\T \Pb)) , (43) 
APsTdHw) = 2 fl« £' A^irolnxXn^l^igaQfc) 



Til 



£Pa &n\ 

(P6|T |n 1 )(Pan 1 |/(A)|gag6) 



£p& — £■ 



ni 



+P (1) (PaPb\l'(A)\QaQb) ({Pa\T \Pa) + (g&|T |g6> - (P&|T |P&)) . (44) 



10 



One can see that the sum of the UV-divergent terms (30), (39)-(44) is zero. 

There are also infrared (IR) divergencies in Eqs. (31)— (38). They arise when the energies of 
the intermediate states n are equal to the energies of the reference states a or b, what leads to the 
factor — iO) 2 . In order to handle these divergencies we introduce a non-zero photon mass fi 
and isolate analytically the terms proportional to In ji. Similar terms arise in Eqs. (45)-(46) and 
should be considered together to yield a finite result. 

C. "Double-vertex" diagrams 

Finally, we consider the diagram D with two vertices inside the self-energy loop, 

i [, (Pan 3 |/H|niQa)(niP6|/(A)|n 2 Q6)(7i2|To|n3) 



SQED 2n J ( £pa _ u _ U£ni ) ( £Qa _ u _ M£n2 ) ( £Qa _ u _ ) 



711,2,3 



The reducible contributions with the third power of u in the denominators are denoted as APgqj^, 

A P SE(7) _ a p SE(/l) a p SE(/2) a p SE(/3) .... 
SQED - A -^SQED + SQED + SQED > ( 4t) ) 

a ^se(/i) = i_ f du ST^ (Pan 2 \I(u) \riiPa) (wjTojwg) 

sqed - 2tt 7 ^ (e Pa - w - M£ ni )(£ Pa - w - «£ n2 ) 

"1,2 

x f + ) (PaPb\I(A)\QaQb) , (47) 

\e Pa —u — ue ni ep a -uj - ue n . 2 J 



{Pan 2 \I(u) \ ni Qa) (mP6|/(A) |n 2 Q6> 



£ Pa -U) - Ue ni )(e Qa -U- UE n2 ] 

x / (Pa|T |Pa) + (Qa|T |Qa) ^ (4g) 



Ep a -u-ue ni e Qa -uj-ue r 

AE Iqe§ = (Pa\Z"(e Pa )\Pa}(Pa\T \Pa}(PaPb\I(A)\QaQb} . (49) 

These contributions are UV-finite. However, they contain IR-divergent terms that should be con- 
sidered together with similar terms from Eqs. (32)-(35). Introducing a non-zero photon mass /i 
and isolating the terms, proportional to In /i we find that the sum of all the IR-divergent terms from 
the terms C, F, H, D and / is finite. 
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D. Numerical evaluation 



Evaluation of the two-electron self-energy correction requires momentum- space calculation 
of the zero- and one-potential terms of the "modified self-energy" contributions (AK?q{^, 
A-^sqed' A-^sqed' an d ^^sqed) an d zero-potential terms of the "modified vertex" contribu- 
tions (AK?q^, AE^q^, and AE^^). For the zero- and one-potential terms of the "modified 
self-energy" contributions we employ the numerical procedure, developed for the self-energy dia- 
gram [51, 52]. 

The "magnetic-vertex" part (with the operator T in the vertex) of the "modified vertex" con- 
tributions is somewhat similar to the vertex part of the one-electron self-energy correction to the 
g factor or the hyperfine splitting. Hence, our treatment of its zero-potential term is based on the 
corresponding calculational procedures developed in Ref. [14] for the g factor and in Refs. [9, 36] 
for the hyperfine splitting. The angular integration, however, required significant generalization, 
in both cases, due to the interelectronic-interaction matrix elements. The "interaction- vertex" part 
(with the interelectronic-interaction operator 1(A) in the vertex) resembles the vertex part of the 
two-electron self-energy correction to the energy levels. Therefore, we have developed the nu- 
merical algorithm for the corresponding zero-potential term on the basis of the one presented in 
Ref. [53]. The main difference is the structure of the angular integrals which is substantially more 
complicated, due to the interaction with the external magnetic field. 

The many-potential terms of the "modified self-energy" and "modified vertex" parts, as well 
as the complete contribution of the "double vertex" part are calculated in the coordinate space. 
Angular integration and summation over intermediate angular momentum projections is carried 
out in the standard way. The summation over the complete spectrum of the Dirac equation at 
fixed angular quantum numbers ^1,2,3 is performed using the dual-kinetic -balance (DKB) approach 
[54] with the basis functions constructed from 5-splines [55]. The infinite summation over ^1,2,3 
is terminated at \k\ = 10-15 and the rest of the sum is evaluated by the least-square inverse- 
polynomial fitting. In order to perform the integration over u we employ two different contours, 
thus performing an additional cross-check. Both of them involve the Wick rotation. The first 
contour is the same as in Ref. [51]. The integration is performed along the imaginary axis. 
Besides, the contributions from the poles arising from the terms with e n < e a ,b must be calculated. 
The advantage of this contour is that the calculation is less time-consuming, in particular, due to 
the fact that Im[J(u;)] = when Re[u] = 0. The second contour was proposed in Ref. [52]. The 
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integration is performed along the line Ch [u — ioo,coo + ioo] and along the half-ellipse Cl, going 
between and u in the lower half -plane. The advantage of this contour is that one does not need 
to investigate the pole-structure of the integrand, which is especially complicated for the diagram 
D. 

In order to check the numerical procedure we have performed the calculation in both Feynman 
and Coulomb gauges for the photon propagator corresponding to the interelectronic interaction. 
The individual terms A-I are presented in both gauges for the g factor of 208 Pb 79+ in Table I 
and for the hyperfine splitting of 209 Bi 83+ in Table EL The data for the g factor demonstrate large 
cancellation of the individual terms, leading to the loss of 2 digits in the total value. Due to this fact 
the uncertainty of our results for the g factor is significantly larger than the one for the hyperfine 
splitting. Moreover, for lower values of Z the convergence of the partial-wave expansion worsen, 
and the resulting accuracy becomes unacceptable. This means that for the g factor in the middle- Z 
region a special treatment of the many potential term is required. 

IV. TWO-ELECTRON VACUUM POLARIZATION 

The diagrams, corresponding to the screened vacuum-polarization correction, are shown in 
Fig. 2. Similar to the case of the screened self-energy we divide the total contribution of these 
diagrams into the irreducible and reducible parts. The reducible part should be considered together 
with the related contributions from the second and the third terms in Eq. (9). The irreducible parts 
of each diagram A - F are denoted by the same letter, AE^^~ F \ while the reducible terms 
are collected into three groups, AE^q^' 11 ' 1 ^. The total correction due to the screened vacuum- 
polarization is given by 

= e £m> p+0 H«r + A^«r + + a«> 

b PQ 

+ ABJSS + AUKS? + A£jS5> + AUggg) . (50) 

Here again P and Q are the permutation operators, interchanging a and b, and the summation over 
b runs over all core electron states. The notation A = SQ b — e Pb will be used below. 

The diagrams A, B and E (so-called 'electric-loop' diagrams) involve the matrix elements of 
the standard Coulomb-field-induced vacuum-polarization potential Uy P . The charge renormaliza- 
tion makes this potential finite, and the expressions for the resulting contributions, Uehling and 
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Wichmann-Kroll potentials, can be found, e.g., in Ref. [56]. The contributions of these diagrams 
are given by 

*%3£ = A « 1) + « 2) , (5D 

AF yP(Ai) _ 9 V ', p nlU oi , >i|T |n 2 )(n 2 Pfe|I(A)|QaQfe) 

SQED - Z 2^ \^ a \ U VP\ n l) 7Z " w » PA> 



ni,2 



Ap vp(A2) _ 9 sr' /Pn{Tre i , >iPfr|/(A)|n 2 Qfr)(n 2 |T |Qa) 

A^sqed -2^(Pa|C/ VP K) (£pa _ £m)(£Qa _ £n2) , (53) 



A^gg ^ 2 E ^^ol^) (na |^| Wl) (^(A)1W) ^ (54) 



SQED - SQED + Ajfc SQED > O^) 

a<<-> = 2 e'(^i^i.».) < T oI " 2>< "'» |j(a)I T > - < 56 > 

„ (£Pa — £ni ) \ £ Pb — £n 2 ) 



Ai ^SQED - 2 2^ V^ a I^VpN) 7- — — 1 • 00 

ft 1,2 

The diagrams of type C (so-called 'magnetic-loop' diagrams) involve the matrix elements of the 
magnetic-field-induced vacuum-polarization potential Uyp. It also requires the charge renormal- 
ization to make it finite. Our present treatment of this potential is restricted to the Uehling (free- 
electron-loop) approximation. The expression for this potential in the Uehling approximation in 
case of the hyperfine interaction can be found in Ref. [57] for the point-dipole nuclear model, and 
in Ref. [47] for the sphere model of the nuclear magnetization distribution. In case of the g factor 
this potential in the Uehling approximation is equal to zero. The contribution of the diagrams of 
type C is given by 

A« = 2 EV«lOh>^^^ • (58) 



£Pa E-m 



The diagram of type F contains the interelectronic-interaction operator modified by the vacuum 
polarization 7 V p (w) • O ur present treatment of this term is restricted to the Uehling approximation. 
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The corresponding expression for Iyp(u) in the Feynman gauge can be found, e.g., in Ref. [58]. 
The contribution of this diagram reads 

A<S = 2 E / (Pa '!° |ril) ^ P& l 7 vp(A)lQaQ&) . (59) 

„ ^Pa &n\ 

The contribution of the diagram D is given by the expression 

A£ s V q P ed = (PaPb\I#(A)\QaQb) , (60) 

where the operator Jyp can be derived within the two-time Green function method. In the present 
work we omit this term, assuming its value is relatively small. 
The reducible contributions G, H and / are given by 

A pVP(G) _ a pVP(Gl) a pVP(G2) a pVPCGS) , fin 
SQED - SQED + SQED + SQED > i 01 -* 



A^ed^ = -2^7 -A{Pa\U^\n 1 ){n 1 \T \Pa){PaPb\I{A)\QaQb) 

+ (Pa|l^ P |Pa)(Pa|T |ni)(n 1 P6|/(A)|QaQ6) 
+ (Pa|To|Pa)(Pa|^ 1 P l«i)(^l^(A)|gag&)} , (62) 



AP s V q P eS 2) = 2 | V' [(Pa|^ 1 pl«i)(^i|7 1 o|Pa)(PaP&|/ / (A)|gag6) 

-(PalU^MimPbll'iA^QaQ^^QblTolQb) - (Pb\T \Pb)^j 

"^v (PaPb \I'(A) \ ni Qb) (wi ITplQa) 



+ ( 



+<Pa|E/&>|Pa) 



^Qa £ni 



, (PaPb\r(A) |Qam> (m |T |Q6) 



(63) 



AP s V qeS 3) = (Pa|C^p|Pa)(PaP6|/"(A)|QaQ6> ((Q6|T |Q6) - (P6|T |P6)J , (64) 
^s V qed = ^ (PaP&|/'(A)|gag&) (<Q6|£/#|Q&) - (P&|£/#|P&>) , (65) 
A ^s V qe2 = \ (PaPb\I YP '(A)\QaQb) ((Qb\T \Qb) - (Pb\T \Pb)^j . (66) 
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The numerical calculations of the screened vacuum-polarization corrections are performed in 
coordinate space employing the finite basis set constructed from the DKB-splines [54]. For the 
electric-loop potential C/ VP , that enters AE^^\ AE^^d, AE$q^, and AE^q^, we employ 
the well-known expression for the Uehling part and the approximate formulas from Ref. [59] for 
the Wichmann-Kroll part. The magnetic-loop potential C/yp (AE$q^ and AE^q^) is taken 
in the Uehling approximation only. In case of the g factor this leads to zero contribution. In 
case of the hyperfine structure we use the expression for the extended nucleus from Ref. [47]. 
For the operator iy P (Ai?sQEE? an ^ A-^sqed) we a l so employ the Uehling approximation, the 
corresponding expressions are taken from Ref. [58]. The calculational procedure has been checked 
utilizing the Feynman and Coulomb gauges for the photon propagator mediating the interelectronic 
interaction. The results in both gauges are presented term by term in Tables I and II for the g factor 
of 208 Pb 79+ and for the hyperfine splitting of 209 Bi 83+ , respectively. 

V. RESULTS AND DISCUSSION 
A. g factor 

The total value of the g factor of a Li-like ion can be written as 

g = g D + Ag int + A^qed + A^qed + Ag nuc . (67) 

Here g D is the Dirac value, Ag int is the interelectronic-interaction correction, Ag QED is the one- 
electron QED correction, A^sqed is the screened QED correction, and Ag nuc is the contribution 
of nuclear effects (finite nuclear size, nuclear recoil and nuclear polarization). To the first order in 
a and l/Z the screened QED correction is defined by the set of the two-electron QED diagrams 
evaluated in the present paper. The operator T in this case reads 

T = ft,, [r x a] • U , (68) 

where /i = |e|/2 is the Bohr magneton, H is the magnetic field directed along the 2-axis. The 
corresponding contribution to the g factor is given by 

A(?sqed = A(?!q ED + A^ ED , (69) 
a #sqed P = a ^sqed P / {noHrrij) , (70) 
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where rrij is the z -projection of the total angular momentum. In Table I the contributions of the 
individual terms to AggQ ED and A#sq ED are presented in the Feynman and Coulomb gauges. The 
contributions of Ai?Jq EB ? and AE^^ are zero in the Uehling approximation that we employed 
for C/yp. The diagram F, which is gauge invariant itself, is calculated only in the Feynman gauge. 
It can be seen from the table that the total results in the different gauges are in a fair agreement with 
each other. Additionally, we estimated the Wichmann-Kroll part of the magnetic-loop diagrams 
taking the one-electron value from Ref. [15], and assuming the same screening ratio as for the 
electric-loop diagrams. 

In Table III the individual terms and the total values of the g factor for Li-like lead 208 pb 79+ and 
uranium 238 U 89+ are presented. The contribution of the screened self-energy correction A(?gQ ED 
amounts to -3.3(2) x 1(T 6 for 208p b 79+ and _ 4 9(9) x 10 -e for 238tj89+ The estmiatec i un- 
certainty of the result is rather large due to the cancellation of the individual terms presented in 
Table I. The values obtained previously with local screening potentials are —3.5(1.2) x 1CT 6 
and —3.1(1.5) x 1CT 6 , respectively [27]. The contribution of the screened vacuum-polarization 
AtfsQED is 1.53(3) x 1(T 6 for 2 ospb 79 + and 2.55(5) x 1(T 6 for 238 U 89+ . The other contributions 
to the g factor presented in Table III were considered in details in our previous studies [26, 27]. 
The accuracy of the g factor value is about 10% better than that from Ref. [27] and is almost 
completely determined by the higher-orders of the interelectronic-interaction correction. 



B. Hyperfine splitting 



The total value of the hyperfine splitting of a Li-like ion can be written as, 

(71) 



AE^l = Ep 



1 1 

A(aZ)(l -S)(l-e) + -B(aZ) + -^C(Z, aZ) + x QED + ^sqed 



where E F is the non-relativistic value of the hyperfine splitting (Fermi energy), A(aZ) is the one- 
electron relativistic factor, 5 and e are the corrections for distributions of the charge and magnetic 
moment over the nucleus, respectively. The interelectronic-interaction corrections of first and 
higher orders in 1/Z are represented by the functions B(aZ) and C{Z,otZ), respectively. The 
term xq E d corresponds to the one-electron QED corrections. The details on these contributions 
can be found in Refs. [40, 47] and references therein. The operator T for the hyperfine splitting 
is given by 

T, = E f gX^. (72) 



iy3 
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The factor G a is defined by the quantum numbers of the valence state, 

_ n 3 (2Z + l)j(j + l) 
° a ~ 2(aZy mj ' (73) 

where n is the principal quantum number, j and rrij are the angular momentum and its projection, 
and / defines the parity of the state. We note that T is the effective one-particle operator, which is 
employed in calculations of various contributions to the hyperfine splitting. The full Hamiltonian 
of the hyperfine interaction is the well-known Fermi-Breit operator. 

To the first order in a and l/Z the screened QED correction ^sqed to the hyperfine splitting is 
given by 

^SQED = ^SQED + ^SQED ? (74) 

SE/VP _ a pSE/VP /p nt -s 
X SQED -AE SQED /E F - ('->) 

In Table IV the screened self-energy and vacuum-polarization corrections to the hyperfine splitting 
are presented for several values of Z in the range Z = 20 - 83. The calculations are performed 
with the Fermi model for the finite nuclear charge distribution. The finite nuclear magnetization 
distribution is introduced via an additional factor F{r) in the operator T [47]. The individual 
contributions to sjIqed an d ^sqed f° r lithiumlike bismuth 209 Bi 80+ are presented in Table II in the 
Feynman and Coulomb gauges. Perfect agreement is found between the total results in the differ- 
ent gauges. This is also true for the other values of Z. We mention, however that the contribution 
of the screened vacuum-polarization diagram F, which is gauge invariant itself, is calculated in the 
Feynman gauge only. For the screened Wichmann-Kroll magnetic-loop part we have employed 
the hydrogenic 2s value from Ref. [57], assuming that it enters with the same screening ratio as 
the Uehling terms. 

In Table V we present the total value of the hyperfine splitting in bismuth in terms of the specific 
difference of the ground state hyperfine splitting in the H-like ion (Is) and in the Li-like ion (2s): 
A'E = AE^ — £Ai^ ls ). It was proposed in Ref. [40] to consider this difference in order to 
overcome the problem of a large uncertainty of the Bohr-Weisskopf (BW) effect, originating from 
the nuclear magnetization distribution. The parameter £ is chosen to cancel the BW correction, 
and the accuracy of the specific difference A'E appears to be much higher than the accuracy of the 
splittings AE 1 - 1 ^ and AE^ themselves. The value of £ = 0.16886 has been found for bismuth, 
taking into account the BW effect on all of the contributions, presented in Table V. For the related 
discussion we refer to our Letter [48]. We only mention here that the rms radius was taken to 

18 



be (r 2 ) 1 / 2 = 5.5211 fm [60], the nuclear spin and parity I T = 9/2—, and the magnetic moment 
H = 4.1106(2)/mt [61]. 

C. Conclusion 

The rigorous evaluation of the screened QED corrections to the g factor and to the hyperfine 
splitting of heavy Li-like ions within ab initio QED approach has been performed. Previously 
developed procedures for the evaluation of the one-electron QED corrections in presence of the 
external magnetic field and of the two-electron QED corrections to the energy levels have been 
associated and generalized. The complete gauge-invariant set of the two-electron self-energy di- 
agrams with external magnetic field has been calculated. The dominant part of the two-electron 
vacuum-polarization correction has been calculated as well. The electric-loop diagrams have been 
evaluated for both Uehling and Wichmann-Kroll parts. The magnetic-loop diagrams have been 
evaluated in the Uehling approximation. These results improve the accuracy of the theoretical 
predictions for the g factor and the hyperfine splitting of heavy ions where stringent tests of the 
bound-state QED effects are feasible. 
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FIG. 1 : Feynman diagrams representing the screened self-energy correction in the presence of an external 
magnetic field. The wavy line indicates the photon propagator and the double line indicates the electron 
propagators in the Coulomb field. The dashed line terminated with the triangle denotes the interaction with 
the magnetic field. 
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FIG. 2: Feynman diagrams representing the screened vacuum-polarization correction in the presence of an 
external magnetic field. The notations are the same as in Fig. 1. 
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FIG. 3: Feynman diagrams representing the two-loop QED corrections in the presence of an external po- 
tential. The notations are the same as in Fig. 1. 
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TABLE I: Contributions of the individual terms to the screened self-energy and vacuum-polarization cor- 
rections to the g factor of Li-like lead 208 pb 79+ . The units are 10 -6 . 





Screened SE 
Feynman Coulomb 




Screened VP 
Feynman Coulomb 


A 


171.7 


173.6 


A 


-34.27 


-34.65 


B 


18.1 


18.5 


B 


0.83 


0.82 


C 


-30.5 


-31.6 


C 






D 


-52.3 


-52.0 


D 






E 


-1.0 


-1.0 


E 


0.16 


0.15 


F 


4.2 


4.2 


F 


-0.01 




G 


-167.3 


-169.2 


G 


34.89 


35.29 


H 


-41.3 


-41.3 


H 






I 


95.3 


95.5 


I 


0.00 




Total SE 


-3.3 


-3.3 


Total(A-I) 


1.60 


1.60 



WK-ml -0.06(3) 
Total VP 1.54(3) 
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TABLE II: Contributions of the individual terms to the screened self-energy and vacuum-polarization cor- 
rections to the hyperfine splitting of Li-like bismuth 209 Bi 80+ in terms of xsqed- 



Screened SE Screened VP 





Pf^vnm an 


cw 1 1 o in h 

\J U l\J 111 u 




T^vnm an 
i t^yiiiiidii 


V- U L11LJ111U 


A 
r\ 




001 


A 


0004881 


0004 8Q 9 


t) 
D 


000^80 


000^ Q 8 


v> 

D 


0009 1 9 8 


0009 1 0^ 


C 


0.001928 


0.001952 


c 


-0.0001691 


-0.0001669 


D 


-0.000936 


-0.000945 


D 






E 


0.000028 


0.000028 


E 


-0.0000031 


-0.0000029 


F 


-0.000174 


-0.000172 


F 


0.0000015 




G 


-0.001298 


-0.001307 


G 


0.0002766 


0.0002749 


H 


0.000331 


0.000331 


H 


0.0000023 


0.0000001 


I 


0.000066 


0.000066 


I 


0.0000000 




Total SE 


0.001109 


0.001109 Total(A-I) 


-0.0005927 


-0.0005927 



WK-ml 0.00005(2) 
Total VP -0.00054(2) 
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TABLE III: Individual contributions to the ground-state g factor of Li-like lead 2U8 pb' 9H ~ and uranium 

238tj89+ 







208 Pb 79+ 


238tj89+ 


Dirac value (point nucleus) 




1.932 002 904 


1.910 722 624(1) 


Finite nuclear size 




0.000 078 58 (13) 


0.000 241 30 (43) 


Interelectronic interaction, ~ 


l/Z 


0.002 148 29 


0.002 509 84 


Interelectronic interaction, ~ 


l/Z 2 and h.o. 


-0.000 007 6 (27) 


-0.000 008 5 (38) 


QED, ~ a 




0.002 411 7(1) 


0.002 446 3 (2) 


QED, ~ c? 




-0.000 003 6 (5) 


-0.000 003 6 (8) 


Screened SE 




-0.000 003 3 (2) 


-0.000 004 9 (9) 


Screened VP 




0.000 001 54 (3) 


0.000 002 55 (5) 


Nuclear recoil 




0.000 000 25 (35) 


0.000 000 28 (69) 


Nuclear polarization 




-0.000 000 04 (2) 


-0.000 000 27 (14) 


Total 




1.936 628 7 (28) 


1.915 905 7 (41) 



TABLE IV: Screened QED corrections to the hyperfine splitting of lithiumlike ions. 



7 ™SE 
^ SQED 


-VP 
X SQED 


^SQED 


20 0.000256 


-0.000116 


0.000140(1) 


30 0.000330 


-0.000131 


0.000199(1) 


40 0.000394 


-0.000155 


0.000238(1) 


50 0.000473 


-0.000186(3) 


0.000287(3) 


60 0.000582 


-0.000241(4) 


0.000340(4) 


70 0.00075 


-0.00033(1) 


0.00042(1) 


83 0.00111 


-0.00054(2) 


0.00057(2) 
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TABLE V: Individual contributions to the specific difference A'E of the hyperfine splittings for bismuth 
209 Bi. The units are meV. 



A£(2s) £AE^ A'E 



Dirac value 844.829 876.638 -31.809 

Interelectronic interaction, ~ 1/Z -29.995 -29.995 

Interelectronic interaction, ~ 1 /Z 2 and h.o. 0.25(4) 0.25(4) 

QED - 5.052 -5.088 0.036 

Screened SE 0.381 0.381 

Screened VP - 0.187(6) - 0.187(6) 

Total -61.32(4) 
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